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The constrained minimal supersymmetric standard model with /i > supplemented by an 'asymp- 
totic' Yukawa coupling quasi-unification condition, which allows an acceptable 6-quark mass, is rein- 
vestigated. Imposing updated constraints from the cold dark matter abundance in the universe, B 
physics, the muon anomalous magnetic moment, and the mass m^ of the lightest neutral CP-even 
Higgs boson, we find that the allowed parameter space is quite limited but not unnaturally small 
with the cold dark matter abundance suppressed only via neutralino-stau coannihilations. The light- 
est neutralino with mass in the range (341 — 677) GeV is possibly detectable in the future direct cold 
dark matter searches via its spin-independent cross section with nucleon. In the allowed parameter 
space of the model, we obtain ruh = (117 — 122.2) GeV. 

PACS numbers: 12.10.Kt, 12.60.Jv, 95.35.+d 



I. INTRODUCTION 

The well-known constrained minimal supersymmetric 
standard model (CMSSM) ilr4], which is a highly pre- 
dictive version of the minimal supersymmetric standard 
model (MSSM) based on universal boundary conditions 
for the soft supersymmetry (SUSY) breaking parame- 
ters, can be further restricted by being embedded in a 
SUSY grand unified theory (GUT) with a gauge group 
containing SU{A)c and SU{2)ji. This can lead [5| to 
'asymptotic' Yukawa unification (YU) [6|, i.e. the exact 
unification of the third generation Yukawa coupling con- 
stants ht, hb, and hr of the top quark, the bottom quark, 
and the tau lepton, respectively, at the SUSY GUT scale 
Mgut- The simplest GUT gauge group which contains 
both S'C/(4)c and SU{2)r is the Pati-Salam (PS) group 
Gps = SU{'i)c X S'C/(2k X SU{2)r 0,9 - for YU within 
50(10), seeRefs. ^,^. 

However, given the experimental values of the top- 
quark and tau-lepton masses (which, combined with YU, 
naturally restrict tan/3 ~ 50), the CMSSM supplemented 
by the assumption of YU yields unacceptable values of 
the 6-quark mass mf, for both signs of the parameter 
II. This is due to the presence of sizable SUSY cor- 
rections [ll| to rub (about 20%), which arise [ll|, [l3 
from sbottom-gluino (mainly) and stop-chargino loops 
and have the sign of /i - with the standard sign conven- 
tion of Ref. [13 ■ The predicted tree-level mb{Mz), which 
turns out to be close to the upper edge of its 95% confi- 
dence level (c.l.) experimental range receives, for /z > 
[^ < 0], large positive [negative] corrections which drive 
it well above [a little below] the allowed range. Conse- 
quently, for both signs of /i, YU leads to an unacceptable 
mb{Mz) with the /i < case being much less disfavored. 
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The usual strategy to resolve this discrepancy is the 
introduction of several kinds of nonuniversalities in the 
scalar 0, [13] and/or gaugino [H, [H sector of MSSM 
with an approximate preservation of YU. On the con- 
trary, in Ref. [ig] , concrete SUSY GUT models based on 
the PS gauge group are constructed which naturally yield 
a moderate deviation from exact YU and, thus, can al- 
low acceptable values of the 6-quark mass for both signs 
of ji within the CMSSM. In particular, the Higgs sector 
of the simplest PS model [3, Q is extended so that the 
electroweak Higgs fields are not exclusively contained in 
a SU{2)l X SU{2)ji bidoublet superfield but receive sub- 
dominant contributions from other representations too. 
As a consequence, a moderate violation of YU is nat- 
urally obtained, which can allow an acceptable 6-quark 
mass even with universal boundary conditions. It is also 
remarkable that the resulting extended SUSY PS models 
support new successful versions [17] of hybrid inflation 
based solely on renormalizable superpotential terms. 

These models provide us with a set of 'asymptotic' 
Yukawa quasi-unification conditions which replace exact 
YU. However, applying one of these conditions in the 
/i < case does not lead [1», Iff] to a viable scheme. 
This is due to the fact that the parameter space allowed 
by the cold dark matter (CDM) requirements turns out 
[ISl Il9| to lie lower than the one allowed by the inclusive 
decay 6 ^- 57 in the ttilsp — A^^ plane, where ttilsp is the 
mass of the lightest sparticle (LSP), which, in our case, is 
the lightest neutralino x and Af^ = (TOf2 — "T'LSp)/"t.lsp 
is the relative mass splitting between the LSP and the 
lightest stau mass eigenstate f2. This result is strength- 
ened by the fact that /Lt < is strongly disfavored by the 
constraint arising from the deviation i5a^ of the measured 
value of the muon anomalous magnetic moment a^ from 
its predicted value af^ in the standard model (SM). In- 
deed, fj, < is defended only at 3a by the calculation 
of af^ based on the r-decay data, whereas there is a 
stronger and stronger tendency at present to prefer the 
e"'"e~-annihilation data for the calculation of a^J^, which 
favor the fj, > regime. Given the above situation, we 



focus here on the /i > case. 

Let us recall that, m this case, the suitable 'asy mptotic' 
Yukawa quasi- unification condition applied [ifl, llal is 



ht : hi, : hr = \1 



II 



|l + 3c| 



(1) 



This relation depends on a single parameter c, which is 
taken, for simplicity, to be real and lying in the range 
< c < 1 . With fixed masses for the fermions of the third 
generation, we can determine the parameters c and tan /3 
so that Eq. (H]) is satisfied. In contrast to the original 
version of the CMSSM ^^, therefore, tan/3 is not a 
free parameter but it can be restricted, within our set- 
up, via Eq. ([T|) to relatively large values. The remaining 
free parameters of our model are the universal soft SUSY 
breaking parameters defined at Mqut, i-e.. 



M- 



1/2, 



rriQ, and Aq, 



(2) 



where the symbols above denote the common gaugino 
mass, scalar mass, and trilinear scalar coupling constant, 
respectively. These parameters can be restricted by em- 
ploying a number of experimental and cosmological re- 
quirements as in Refs. [la, [l3|- In view of the expected 
data from the Large Hadron Collider (LHC), it would be 
worth to retest our model against observations using the 
most up-to-date version of the available constraints. 

We exhibit the cosmological and phenomenological re- 
quirements which we considered in our investigation in 
Sec, mi and we restrict the parameter space of our model 
in Sec. IIIII Finally, we test our model from the perspec- 
tive of the CDM direct detection experiments in Sec. IIVI 
and summarize our conclusions in Sec. |Vl 



II. COSMOLOGICAL AND 
PHENOMENOLOGICAL CONSTRAINTS 

In our investigation, we integrate the two-loop renor- 
malization group equations for the gauge and Yukawa 
coupling constants and the one-loop ones for the soft 
SUSY breaking parameters between A/gut and a com- 
mon SUSY threshold Msusy — (™ti"i(J^/^ (ii^2 are 
the stop mass eigenstates) determined in consistency 
with the SUSY spectrum. At Msusy, we impose ra- 
diative electroweak symmetry breaking, evaluate the 
SUSY spectrum employing the publicly available calcula- 
tor SOFTSUSY HO], and incorporate the SUSY corrections 
to the b and r mass [l2] . The corrections to the r-lepton 
mass rur (almost 4%) lead [l^ [31 to a small decrease 
of tan/3. From Msusy to Mz, the running of gauge and 
Yukawa coupling constants is continued using the SM 
renormalization group equations. 

The parameter space of our model can be restricted 
by using a number of phenomenological and cosmologi- 
cal constraints. We calculate them using the latest ver- 
sion of the publicly available code micrOMEGAs [21|. We 
now briefiy discuss these requirements - for similar re- 
cent analyses, see Ref. [H for CMSSM or Refs. [11 [H 
for MSSM with YU. 



a. SM Fermion Masses. The masses of the fermions 
of the third generation play a crucial role in the determi- 
nation of the evolution of the Yukawa coupling constants. 
For the &-quark mass, we adopt as an input parameter in 
our analysis the MS 5-quark mass, which at la is [24] 



/ \MS , in 



0.06 ^'^^ ■ 



(3) 



This range is evolved up to Mz using the central value 
as{Mz) = 0.1184 [2J] of the strong fine structure con- 
stant at Mz and then converted to the DR scheme in 
accordance with the analysis of Ref. |25j . We obtain, at 
95% C.I., 



2.745 < TO6(Mz)/GeV < 3.13 



(4) 



with the central value being mb{Mz) = 2.84 GcV. For 
the top-quark mass, we use the central pole mass (Mt) 
as an input parameter [26| : 



Mt = 173 GeV 



mtimt) = 164.6 GeV (5) 



with mt{mt) being the running mass of the t quark. We 
also take the central value mT{Mz) = 1.748 GeV ^] of 
the DR tau-lepton mass at Mz- 

b. Cold Dark Matter Considerations. According to 
the WMAP results ^, the 95% c.l. range for the CDM 
abundance is 



rtcDuh^ = 0.1126 ±0.0072. 



(6) 



In the context of the CMSSM, the LSP can be the lightest 
neutralino x and naturally arises as a CDM candidate. 
We require its relic abundance JIlsp^^ in the universe 
not to exceed the upper bound derived from Eq. ([6]) - 
the lower bound is not considered since other production 
mechanisms 28| of LSPs may be present too and/or other 
CDM candidates [23, 30,] may also contribute to flcDMh"^- 
So, at 95% C.L, we take 



^LSP^^ ^ 0.12. 



(7) 



An upper bound on ttilsp (or m^) can be derived from 
Eq. ([7]) since, in general, JIlsf^^ increases with tolsp- 
The calculation of JIlsp^^ in micrOMEGAs includes accu- 
rately thermally averaged exact tree-level cross sections 
of all the possible (co)annihilation processes [3|,|3l|, treats 
poles [j, [la, 133 properly, and uses one-loop QCD and 
SUSY QCD correctfons [ll [H HI Hi to the Higgs de- 
cay widths and couplings to fermions. 

c. The Branching Ratio BR {b — >■ s"f) ofb-^ sj. The 
most recent experimental world average for BR(6 — ?> 57) 
is known fsBl] to be (3.52 ± 0.23 ± 0.09) x 10"^ and its 
updated SM prediction is (3.15 ± 0.23) x lO"'' ^. Com- 
bining in quadrature the experimental and theoretical er- 
rors involved, we obtain the following constraints on this 
branching ratio at 95% c.l.: 



2.84 X 10"* < BR( 



57) < 4.2 X 10" 



(8) 



The computation of BR(6— >-S7) in the micrOMEGAs 
package presented in Rcf. [SJ] includes [33l next-to- 
leading order (NLO) QCD corrections to the charged 
Higgs boson (H^) contribution, the tan/3 enhanced con- 
tributions, and resumnied NLO SUSY QCD corrections. 
The H contribution interferes constructively with the 
SM contribution, whereas the SUSY contribution inter- 
feres destructively with the other two contributions for 
^ > 0. The SM plus the H"^ and SUSY contributions 
initially increases with ttilsp and yields a lower bound 
on TTiLSP from the lower bound in Eq. ([S]). (For higher 
values of tolsp, it starts mildly decreasing.) 

d. The Branching Ratio BR (i?s — >■ fi'^fi") of Bg — >■ 
fi^ fi^ . The rare decay B^ — )■ fJ-^ IJ-^ occurs via Z pen- 
guin and box diagrams in the SM and, thus, its branch- 
ing ratio is highly suppressed. The SUSY contribution, 
though, originating [3^, [33] from neutral Higgs bosons in 
chargino-, H -, and VF -mediated penguins behaves as 
tan^ f3/m\ and hence is particularly important for large 
tan/3's. We impose the following 95% c.l. upper bound: 



BR B, 



fJ' 



>") < 5.8 X 10" 



(9) 



as reported [40| by the CDF collaboration. This bound 
implies a lower bound on ttilsp since BR(i?s — s> ^i^fi^) 
decreases as ttilsf increases. 

e. The Branching Ratio BR (i?„ -^ tv) of i?„ — s> tv. 
The purely leptonic decay i?„ -^ tv proceeds via W^- 
and JT^-mediated annihilation processes. The SUSY 
contribution, contrary to the SM one, is not helicity sup- 
pressed and depends on the mass m//± of the charged 
Higgs boson since it behaves [33,|41j as tan'* /3/m'^± . The 
ratio R (B„ -^ tv) of the CMSSM to the SM branching 
ratio of B^ — > tv increases with ttilsp and approaches 
unity. It is to be consistent with the following 95% c.l. 
range |35| : 



0.52 <R(B„ -^Tv) < 2.04 



(10) 



A lower bound on ttzlsf can be derived from the lower 
bound in this inequality. 

/. Muon Anomalous Magnetic Moment. The quan- 
tity (5ap, which is defined in Sec. U can be attributed 
to SUSY contributions arising from chargino-sneutrino 
and neutralino-smuon loops. The relevant calculation is 
based on the formulas of Ref. |42| . The absolute value of 



the result decreases as ttilsp increases and its sign is pos- 
itive for /u > 0. On the other hand, the calculation of a™ 
is not yet stabilized mainly because of the ambiguities in 
the calculation of the hadronic vacuum-polarization con- 
tribution. According to the most up-to-date evaluation of 
this contribution in Ref. [ijl, there is still a discrepancy 
between the findings based on the e"'"e~-annihilation data 
and the ones based on the r-decay data. Taking into ac- 
count the more reliable calculation based on the e+e^ 
data and the experimental measurements [4J] of a^ , we 
obtain the following 95% c.l. range: 



12.7 X 10"^° < (5a^ < 44.7 x 10" 



^Lsph =0-12 



X 



BR(b^SY)< 2.84x10 
BR(B^-*|jV)> 5.8x10"^ 
R(B^^Tv)<0.52 

5a > 12.7x10"'° 

V 

5a =2.9x10"'° 

v 

m, < 114.4 GeV 

h 

nLgph'<0.12 & 
5a > 2.9x10"'° 



FIG. 1: Summary of the conventions adopted in Figs. [2] and 
[3] for the various restrictions on the parameters of the model. 



The T-decay based calculation, on the other hand, yields 
the following 95% c.l. range: 



2.9 X 10"*° < Son < 36.1 X 10~*°. 



(12) 



A lower [upper] bound on mLSP can be derived from the 
upper [lower] bound in Eqs. ([TT|) and ([T2|). As it turns 
out, only the upper bound on ttilsp is relevant in our 
case. Taking into account the aforementioned computa- 
tional instabilities, we will impose the less stringent up- 
per bound on tolsp from the r-decay based calculation. 
However, we will also depict the more stringent bound 
from the e~'"e~-annihilation data for comparison. 

g. Collider Bounds. For our analysis, the only rele- 
vant collider bound is the 95% c.l. LEP bound [45] on 
the lightest CP-even neutral Higgs boson mass 



m,, 



, > 114.4 GeV, 



(13) 



10 



(11) 



which gives a lower bound on tolsp- The calculation of 
ruh in the package SOFTSUSY [lo] includes the full one- 
loop SUSY corrections and some zero-momentum two- 
loop corrections [4^. The results are well tested [43] 
against other spectrum calculators. 



III. RESTRICTIONS ON THE SUSY 
PARAMETERS 

Imposing the requirements described above, we can de- 
lineate the allowed parameter space of our model. The 
predicted mass spectra are possibly relevant for the LHC 
searches. Throughout our investigation, we consider 
the central values for the SM parameters Mj, mi,{Mz), 
niriMz), and as{Mz). We adopt the following conven- 
tions for the various lines and regions in the relevant fig- 
ures (Figs. [2] and [3]) - see Fig.[T] 

• on the solid black line, Eq. ([7]) is saturated, 
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FIG. 2: Restrictions in the Afi/2 — mo plane for various values of Ao/Afi/2 indicated in the graphs. The conventions adopted 
are described in Fig. [1] 



• the light gray region is cosmologically excluded 
since it predicts charged LSP, 

• the dark gray region is excluded by the lower bound 
inEq. ©, 

• the gray region is excluded by Eq. (j9]), 

• the yellow region is excluded by the lower bound in 

Eq. (nnD, 

• the vertically hatched region is favored by the lower 
bound in Eq. (fTTI) . 

• on the dotted black line, the lower bound in 
Eq. p^ is saturated, 

• the red region is excluded by Eq. P^ . 

• the horizontally hatched region is allowed by both 
Eq. ([7]) and the lower bound in Eq. ([T^ . 

Note that the upper bounds in Eqs. ©, dTU]), (HH), and 
(IT^ do not restrict the parameters of our model. 



We present the restrictions from all the requirements 
imposed in the M1/2 — "tuq plane for Ao/Afi/2 = 0, 1, —1, 
and —2 in Fig. [2j We remark that the lower bound on 
Afi/2 comes from Eq. ^ for ylo/A/1/2 = 0, -1, and -2 
and from the lower bound in Eq. p^ for AqJMxi^ = 1- 
Also, from the relevant data, we observe that the lower 
bound in Eq. ((TO)) is fulfilled for the mass of the CP- 
odd Higgs boson tua — 520 GeV and almost indepen- 
dently of the other parameters. Finally, note that, for 
Ao/Mi/2 = —1 and —2, the bound in Eq. p^ is violated 
for Ml/2 < 400 GeV and, consequently, does not appear 
in the relevant diagrams. 

The constraint in Eq. ([7]) is, in general, satisfied in two 
well-defined distinct regions in the diagrams of Fig. [21 In 
particular, 

• the region to the left of the almost vertical part 
of the line corresponding to the upper bound on 
Ml/2 from Eq. ([7]), where the LSP annihilation via 
the s-channel exchange of a CP-odd Higgs boson 
A is by far the dominant (co) annihilation process. 
However, this region is excluded by the constraints 



in Eqs. ([9]) and (TTOl) . On the other hand, it is weh 
known - see e.g. Refs. [l,[ia| - that this region is ex- 
tremely sensitive to variations of mb{Mz). Indeed, 
we find that as inb{Mz) decreases, the A-boson 
mass ruA increases and approaches 27tilsp. The A- 
pole neutralino annihilation is then enhanced and 
f^LSP^^ is drastically reduced causing an increase 
of the upper bound on M^ /2. However, even if we 
reduce mb{Mz), we do not find any ^-pole annihi- 
lation region which is allowed by the requirements 
of Eqs. dSJ and (Uni). 



• the narrow region which lies just above the light 
gray area with charged LSP, where bino-stau coan- 
nihilations [3, l3l| take over leading to a very pro- 
nounced reduction of Vl^sph^- A large portion of 
this region survives after the application of the re- 
quirements in Eqs. ([9]) and ([T0|) and constitutes the 
overall allowed parameter range of our model for 
the given Aq. To get a better understanding of this 
area, we can replace the parameter ttiq by the rel- 
ative mass splitting Af^ between the LSP and the 
lightest stau, defined in SeclH We observe that the 
overall allowed region requires Af^ < 0.025. It is 
evident from Fig. [2] that the slope of the bound- 
ary line with Af^ = increases as A0/M1/2 moves 
away from zero in both directions. Note that this 
slope in our model turns out to be larger than the 
one obtained in other versions of the CMSSM ~ cf. 
Ref. [3] - with lower values of tan/3. As a conse- 
quence, small variations of toq or A/x/2 lead, in our 
model, to more drastic variations in Af^. 

Finally, we note that the more stringent upper bound 
on AIi/2 from the lower bound in Eq. pTJ) is not sat- 
isfied for the values taken for Aq/Mi/2 in Fig. [21 with 
the values Aq/Mi/2 = and 1 being much more favored. 
On the other hand, the lower bound in Eq. P^ is ful- 
filled in the whole allowed region for Aq/Mi/2 = ^nd 1 
whereas, for Aq/AIi/2 — —1 and —2, it imposes an upper 
bound on M1/2 which overshadows the bound on M1/2 
from Eq. ^ . Since the saturation of the lower bound in 
Eq. (fT2)) occurs for Af^ < 0.01, the portion of the dotted 
black line - see Fig. [1] - which connects the black solid 
line with the boundary of the gray area is not visible in 
the relevant panels of Fig. [5] 

To get a better idea of the allowed parameter space, 
we focus on the coannihilation regime and construct the 
allowed region in the M1/2 — Aq/Mi/2 plane. This is 
shown in Fig. [3J where we depict the restrictions on the 
parameters from the various constraints for Af^ — 0. 
This choice ensures the maximal possible reduction of 
^LSP^^ due to the x~^2 coannihilation. So, for A^^ = 0, 
we find the maximal M1/2 or tolsp allowed by Eq. ([7]) 
for a given value of Ao/Mi^2- We observe that, for 
-0.8 < Ao/Mi/2 ^ 3 [-2.55 < A0/M1/2 < -0.8 and 
3 ^ Aq/AIi/2 ^ 3.21] the overall upper bound on Afi/2 or 
"T-LSP is derived from the bound in Eq. ([T]) [lower bound 
in Eq. ^^]. We find that, for A0/M1/2 < 0, processes 
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FIG. 3: Restrictions in the M1/2-A0/M1/2 plane for A^j = 
following the conventions of Fig. [U but with the horizontally 
hatched region not extended to areas excluded by other con- 
straints. 



with f2f2* in the initial state and W^W^ , W^H^ in 
the final one become more efficient (with a total contri- 
bution to the effective cross section of about 14 to 21% 
as AqIMxI2 decreases from to -2.55) and so coanni- 
hilation is strengthened and tolsp's larger than in the 
Ao/Mi/2 > case are allowed. The overall maximal 
Ml/2 - 1495.4 GeV or mLSP ^ 677 GeV is encoun- 
tered at ylo/Afi/2 — —0.8. On the other hand, for 
-2.55 < Ao/Mi/2 ^ 0.7 [0.7 < Ao/A/1/2 < 3.21] the 
lower bound on M1/2 or r?T.LSP is derived from the bound 
in Eq. ^ [lower bound in Eq. ([10])]. The overaU al- 
lowed lowest Afi/2 '^^ 11\.22 GeV or mLSP ^ 341 GeV 
is encountered at ^o/A^i/2 — 0.7. Let us remark that 
the more stringent upper bound on Mxi2 from the lower 
bound in Eq. (fTTj) is not satisfied in the allowed region 
of our model, since there is no common region between 
the horizontally and the vertically hatched areas for any 
Ao/A/1/2. However, for < ^o/-^^i/2 ^ 1; these ar- 
eas are quite close to each other. Note that increasing 
Af2 within its allowed range — 0.025 does not alter the 
boundaries of the various constraints in any essential way, 
except the solid line which is displaced to the left so that 
the allowed area shrinks considerably. 

The deviation from YU can be estimated by defining 
p^ the relative splittings 8hb and Sh^ at Mqut through 
the relations: 



5hh 



2c 



= — (5ft-r = 



ht 



(14) 



In the allowed (horizontally hatched) area of Fig. [31 the 
ranges of the parameters c, J/i,-, ^^ib, and tan/3 are 



0.149 <c< 0.168, 

0.26 < 8hr == -bhb < 0.29, 

56.3 <tan^< 57.7. 



(15) 




TABLE I: Input and output parameters, masses of the spar- 
ticles and Higgs bosons, and values of the low energy observ- 
ables of our model for four values of Ao/Mi/2- Recall that 
1 pb~2.6 X 10-" GeV-^ 
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FIG. 4: The allowed (horizontally hatched) region in the 
tolsp — ruh plane for Afa — 0. We also depict the curves 
corresponding to various values of A0/M1/2, indicated on 
them. The dark points on the boundary correspond to 
Ao/Mi/2 = -2.55, -0.8, 3, 3.21, and 0.8 starting from the 
point at the top of the allowed area and moving clockwise. 



Let us underline that, although the required deviation 
from YU is not so small, the restrictions from YU are 
not completely lost since tan/3 remains large - close to 
60 - and the deviation from YU is generated in a GUT- 
inspired well-motivated way. 

Taking into account the results depicted in Fig. [31 we 
can make predictions for the sparticle and the Higgs bo- 
son spectrum of our model, which may be observable at 
the LHC. In Table U we list the model input and out- 
put parameters, the masses in GeV of the sparticles - 
neutralinos x, X21 Xs: X%i charginos xi^, xj, gluinos g, 
squarks ti, ^2, ^1, ^2, ul, ur, d^, d^, and sleptons fi, 
T2, Vr, cl, en, i>e ~ and the Higgs bosons {h, H, H^, A), 
and the values of the various low energy observables for 
Ao/Mi/2 = 0, ±1, and —2 and for the lowest possible 
Ml/2 in each case adjusting Af^ so as J^lsp^^ — 0.11. 
Note that we consider the squarks and sleptons of the 
two first generations as degenerate. From the values of 
the various observable quantities, it is easy to verify that 
all the relevant constraints are met. We also included in 
Table U predictions for the possible direct detection of the 
LSP using central values for the hadronic inputs f^ or 



AP-see SeclWl 

For the lowest masses of the Higgs and sparticle spec- 
trum (jrih and ?nLSp)i 'we present even more explicit pre- 
dictions in Fig.m where we depict the allowed m-^'s ver- 
sus mLSP for Af2 — and Aq/Mi/2 = 0, ±1, and ±2. As 
can be seen from Fig. [31 the lower limits on the solid lines 
for Ao/Mi/2 = 0, -1, and -2 [A0/M1/2 = 1 and 2] are 
found from the bound in Eq. © [lower bound in Eq. ([TU])] 
- see also Table I. On the other hand, the upper limits of 
the sohd lines for A0/M1/2 = 0, 1, and 2 [A0/M1/2 = -1 
and —2] are found from the bound in Eq. ((7|) [lower bound 
in Eq. (H^j)]. The approximate overall allowed area in the 



Input parameters 


Ao/Mi/2 





1 


-1 


-2 


c 


0.161 


156 


0.165 


0.168 


Ml/2/Gev 


825.7 


776.06 


927.25 


1041.8 


mo/GcV 


665.4 


687.5 


943.1 


1466.8 


Output parameters 


tan/3 


57 


56.8 


57.4 


57.7 


ht(MGUT) 


0.58 


0.58 


0.58 


0.58 


1005/1, (Mgut) 


27.7 


26.9 


28.3 


28.7 


/x/GeV 


925.8 


804 


1170 


1505 


Af,(%) 


2.46 


2.45 


2.13 


1.52 


Masses in GeV of sparticles and Higgs bosons 


X 


365.7 


342.7 


413.2 


467.5 


X§ 


705 


656 


802 


909 


X"s 


927 


807 


1170 


1502 


xl 


940 


827 


1177 


1506 


xf 


940 


827 


1177 


1506 


xf 


705 


656 


802 


909 


9 


1916 


1813 


2145 


2412 


ii 


1585 


1530 


1752 


1980 


h 


1383 


1352 


1506 


1666 


61 


1578 


1526 


1752 


2008 


62 


1498 


1454 


1670 


1916 


UL 


1052 


1762 


2134 


2585 


Ur 


1011 


1694 


2054 


2503 


dL 


1055 


1764 


2135 


2586 


dR 


1006 


1764 


2045 


2494 


fi 


777 


754 


956 


1283 


f2 


374.7 


351.1 


422 


474.6 


Vt 


756 


738 


939 


1272 


e.L 


880 


875 


1142 


1635 


e-R 


740 


752 


1010 


1523 


Ve 


876 


871 


1139 


1633 


h 


118.1 


117 


119.7 


121.3 


H 


584 


519 


668 


747.6 


H^ 


591 


527 


674 


752.8 


A 


585 


520 


669 


748 


Low 


energy observables 




10^BR(&->S7) 


3.32 


3.41 


3.32 


3.37 


10**BR {Bs -> ^J^^i^ 


') 5.76 


5.3 


5.78 


5.8 


R {Bu -^ Tv) 


0.61 


0.52 


0.69 


0.74 


W^°5a^, 


10.6 


11.6 


6.9 


3.9 


JIlsp/i 


0.11 


0.11 


0.11 


0.11 


a|;,/10-9pb 


0.536 


1.1 


0.2 


0.076 


4?/io-^pb 


1.96 


4.1 


0.6 


0.2 



tolsp — TT^h plane is hatched. Shown are also the bound- 
ary points of this region at Aq/AIi/2 — —2.55, —0.8, 3, 
3.21, and 0.7 starting from the point at the top of the 
allowed area and moving clockwise. 

As one can see from Fig. |4l mh increases with ttilsp 
and as ^o decreases. Since the maximum allowed ttilsp 
from the bound in Eq. ((7]) or the lower bound in Eq. (fT2|) 
is achieved at Af^ ~ for given Aq, we conclude that 
the maximum possible allowed rrih can be obtained for 
Af2 ~ 0. On the other hand, the minimum possible al- 
lowed TTifi practically coincides with its value for Af^ ~ 
since variation of Af^ within the values allowed by Eq. ([7]) 
causes minor modifications of mh for fixed Mi/2. For 
Ao = 0, we find 826.4 < Mi/2/GeV < 1348.9 or 365.9 < 
TOLSp/GeV < 607.4 and 118.1 < mh/GeY < 120.6. The 
overall minimum [maximum] rrih is 117.03 [122.2] ob- 
tained at Ao/Mi/2 ~ 1 [Ao/Mi/2 ~ -2.55] for M1/2 = 
776.6 GeV [M1/2 = 1106.6 GeV] or mLSP ^ 343.1 GeV 
[mLSP -498.3 GeV]. 



IV. CDM DIRECT DETECTION 

As we have shown, our model possesses a limited and 
well-defined range of parameters allowed by all the rele- 
vant cosmological and phenomenological constraints. It 
would be, thus, interesting to investigate whether the 
predicted LSPs in the universe could be detected in the 
current or planned direct CDM searches |48l45(j ]. which 
look for evidence of weakly-interacting massive particles 
through scattering on nuclei. The quantities which are 
conventionally used in the recent literature for comparing 
experimental results and theoretical predictions are the 
spin-independent (SI) and spin-dependent (SD) lightest 
neutralino-proton (x — p) scattering cross sections cr|p 
and cr^, respectively. 

These quantities are calculated by employing the rel- 
evant routine of the micrOMEGAs package [5l| based on 
the full one-loop treatment of Ref . 52] , which happens to 
agree with the tree-level approximation [5j| for the val- 
ues of the SUSY parameters encountered in our model. 
Following the approach of Refs. [5l|,[53|, we calculate the 
scalar form factors for light quarks in the proton /^ 

(with q = u,d,s), needed for the calculation of cr^l,, via 
the formulas: 



JTd " 



2a, 



N 



1 



n-p _^ ^^U J^U pp 

fP _ VOttN 

JTs — 



(1 + ^) 



rud 



(16a) 

(16b) 
(16c) 



Here we take for the mass of the proton nip = 0.939 GeV 
and for the light quark mass ratios 



0.553 ±0.043 and 



rrid 



rrid 



18.9 ±0.8, 



(17) 



whereas the ratio Bu/Bd is evaluated from 

Bu _ 2z-{z- l)y 
Bd^ 2 + (z - l)y 



(18) 



with z = 1.49. The uncertainties in z and the quark mass 
ratios are negligible compared to the uncertainties in the 
pion-nucleon sigma term a^^isi and the fractional strange 
quark content of the nucleon y, for which recent lattice 
simulations suggest ^] that, at 68% c.L, 



CTttAt = 53 



+21.1 
-7.3 



MeV and y = 0.030 



+0.017 
0.018- 



(19) 



Taking into account the relations above, we find the fol- 
lowing 1(7 ranges for the /|^„'s: 



fP 

JTu 

fP 

fP . 
JTs 



024+°°™^ 

U.UZ^_Q QQ32, 

029+°°i2 

U.UZ»_g QQ42, 



0.021 



-1-0.025 
-0.013- 



(20a) 
(20b) 
(20c) 



Note that f^^ turns out to be considerably smaller than 
its older value - cf. Ref. [l^ - reducing thereby the 
extracted afl,. 



XP' 



For the calculation of o'lJ^', the relevant axial- vector 



XP 



form factors for light quarks in the proton A^ (with q 
u, d, s) are taken to lie in their Icr ranges [55j: 



A^ = 


+0.842 ±0.012, 


(21a) 


AS- 


-0.427 ±0.013, 


(21b) 


A?- 


-0.085 ±0.018. 


(21c) 



Taking the central value of JIlsp^^ in Eq. ©, but al- 
lowing the hadronic inputs /^ or A^ to vary within their 
ranges in Eqs. (I20ap - (|20cp or (|21ap - (|21cp . respectively, we 
derive the dark gray, gray, and light gray hatched bands 
in the tolsp — cr|p or rriLSP — c|p plane corresponding to 
Ao/Mi/2 = 0, 0.7, and -0.8, respectively - Fig. [5l The 
selected values of Aq/Mi^2 allow us to cover the whole 
range of the allowed tolsp's in our model - cf. Fig. H) 
The bold solid lines in the middle of the bands of the 
left panel in Fig. [5] correspond to the central values of 
the /x„'s. We used the central value of i^LSP^^ since, 
as it turns out, for fixed mLSP, cr^p and a^ are almost 
insensitive to the variation of J^lsp^^ within the range 
of Eq. ^ - or, equivalently, to the required variation of 
Af2 . The width of the bands is almost exclusively due 
to the variation of /^ or A^. As a consequence, the 
bands in the tolsp — cr-J, plane are wider than those in 
the mLSP — cr|p plane due to the larger uncertainties in- 
volved in the determination of the /x„'s. In the left panel 
of Fig. [5l we depict by a dashed green line the recently 
announced 49] upper bound on cr|p from XENON which 
is slightly lower than the one from CDMSII [48|, which 
is not included in the panel. We also draw with dot- 
ted red lines the projected sensitivities of SuperCDMS 
at Soudan and SNOLAB [5g] - from top to bottom. Our 
model can be ultimately tested by XENON- 1 ton, whose 
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FIG. 5: The SI and SD x~ P cross sections (t|p and o-^^ , respectively, versus ttilsp for various Ao/Mi/2's indicated in tiie 
grapiis. The bold solid lines in the left panel are derived by fixing Qi^sph^ and /^ to their central values in Eqs. (|6]) and 
(|20a|l - (|20c|l . whereas the hatched bands in both panels by allowing the hadronic inputs f^ or A^ to vary in their ranges in 
Eqs. (|20a|l - (|2Uc|l or (|21ap - (|21cp . The present and planned sensitivity limits of the various experimental projects are also depicted 
by dashed and dotted lines, respectively. 



planned sensitivity [56|, depicted by a green dotted line, 
covers almost the whole available parameter space of the 
model. On the contrary, as can be easily deduced from 
the right panel of Fig. [SJ a^ in our model lies well below 
the sensitivity of IceCube [501 (assuming neutralino an- 
nihilation into VK+VF~) - depicted by a dashed red line 

- and the expected limit from the large DMTPC detec- 
tor [56|], denoted by a dotted green line. Therefore, the 
LSPs predicted by our model can be detectable in the 
future projects which will release data on cr|p. Further- 
more, the overall upper bound on tolsp found in Sec. IIIII 

- "T-LSP ^ 677 GeV - implies lower bounds on a^ and 
ct|^. Namely, 

afp > 4.3 (3.6) X 10"" pb and a^^ > 1.5 (1.4) x 10"^ pb, 

(22) 
where the bounds in parentheses are derived by allowing 
the /xo's and A^'s to vary within la. Needless to say 
that the low values of a^l, and cr|p obtained here are due 
to the fact that we use universal 'asymptotic' gaugino 
masses and, thus, the LSP is an almost pure bino, as in 
every version of the CMSSM. 



CONCLUSIONS 



the CDM abundance in the universe, B physics (6 -^ 57, 
Bs -^ ^J■^fJ■~, and i?„ — ?► tu), (5a^, and mh- We showed 
that our model possesses a limited but not unnaturally 
small range of parameters which is consistent with all 
these requirements. Namely, the constraint arising from 
CDM considerations can be satisfied simultaneously with 
all the other constraints thanks to the drastic reduction of 
the LSP relic density by neutralino-stau coannihilations. 
For Aq = 0, we find 365.9 < mLSp/GeV < 607.4 and 
118.1 < TO/i/GeV ^ 120.6, whereas, in the overall allowed 
region of our model, we have —2.55 < A0/M1/2 ^ 3.21 
with 341 < mLSp/GeV < 677 and 117 < nih/GeV < 
122.2. Almost all the allowed parameter space of our 
model will be accessible in future CDM direct experi- 
ments which look for SI cross sections between neutralino 
and proton. 

It is worth mentioning that the present investiga- 
tion constitutes an improved version of the analysis in 
Ref. [Ig. The consideration of the constraints from 
BR {Bg -^ M^M ) and R {B^ — > ri^), the updated experi- 
mental results for all the other constraints, and the eval- 
uation of the particle spectrum employing SDFTSUSY are 
the main improvements in this work. The results ob- 
tained are significantly different from the previous ones. 



We performed a revised scan of the parameter space 
of the CMSSM with /i > applying a suitable Yukawa 
quasi-unification condition predicted by the SUSY GUT 
model of Ref. [la] , which has been constructed in order 
to remedy the &-quark mass problem arising from exact 
Yukawa unification and universal boundary conditions. 
We took into account updated constraints from collider 
and cosmological data. These constraints originate from 



Note Added 

While this work was under completion, we became 
aware of Ref. [53] , where the CMSSM with Yukawa quasi- 
unification is also analyzed. Although our results as re- 
gards tan/3, c, cr|p, and a^ are similar, there are large 
discrepancies as regards the CMSSM mass parameters 
and, consequently, the mass spectrum. In particular, 



the ratio {niA — 2?71lsp)/2tolsPj which determines the 
strength of the A-pole effect in reducing r^LSP^^i is not 
allowed to be lower than 0.2 in our case and, thus, this 
effect is excluded. Indeed, the portion of the parameter 
space allowed by Eq. (O due to A-pole neutralino anni- 
hilations is excluded by the B-physics constraints in our 
analysis - contrary to the findings of Ref. [53|. These 
discrepancies can be possibly attributed to the fact that 
we use different numerical routines for the calculation of 
both the SUSY spectra and the low energy observables. 
It is well known [58| that the predictions of the various 
SUSY spectrum calculators do not coincide in the large 
tan/3 regime. Our results as regards the implementation 
of the electroweak symmetry breaking are consistent with 
our initial investigation in Ref. |16l | . 
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